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Abstract 



A search for squarks in i?-parity violating supersymmetry is performed in e'^p collisions 
at HERA at a centre of mass energy of 300 GeV, using HI data corresponding to an inte- 
grated luminosity of 37 pb~^. The direct production of single squarks of any generation in 
positron-quark fusion via a Yukawa coupling A' is considered, taking into account i2-parity 
violating and conserving decays of the squarks. No significant deviation from the Standard 
Model expectation is found. The results are interpreted in terms of constraints within the 
Minimal Supersymmetric Standard Model (MSSM), the constrained MSSM and the mini- 
mal Supergravity model, and their sensitivity to the model parameters is studied in detail. 
For a Yukawa coupling of electromagnetic strength, squark masses below 260 GeV are ex- 
cluded at 95% confidence level in a large part of the parameter space. For a 100 times 
smaller coupling strength masses up to 182 GeV are excluded. 
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1 Introduction 



Supersymmetry (SUSY) relates elementary fermions and bosons and protects the mass of the 
Higgs boson from acquiring unnaturally large radiative corrections. SUSY is often considered 
an ingredient of a fundamental theory beyond the Standard Model (SM). It is thus actively 
searched for in current experiments. 

The ep collider HERA which provides both baryonic and leptonic quantum numbers in the 
initial state is ideally suited to search for new particles possessing couplings to an electron-quark 
pair. Such particles could be squarks, the scalar SUSY partners of quarks, in models where R- 
parity, a discrete symmetry related to lepton and baryon number conservation, is violated (^p). 
These squarks could thus be resonantly produced at HERA via the fusion of the initial state 
positron of energy 27.5 GeV with a quark coming from the incident proton of energy 820 GeV, 
up to the centre of mass energy ~ 300 GeV. 

In this paper, a search is performed for squarks that are singly produced via an coupling, 
considering both decays and decays via gauge couplings involving neutralinos, charginos or 
gluinos. The data were taken from 1994 to 1997 and correspond to an integrated luminosity of 
37 pb^^. This analysis extends the searches for eg resonances previously performed by HI [jT]] 
using the same data sample by considering specific squark decay modes, and supersedes earlier 
published dedicated squark searches [|^, |^] which were based on ~ 13 times less data. 



2 Phenomenology 

The most general SUSY theory which preserves the gauge invariance of the Standard Model 
allows for Yukawa couplings between two known SM fermions and the scalar SUSY partner 
of a quark (a squark q) or of a lepton (a slepton /). Such couplings induce violation of the 
i?-parity defined as Rp = (— i)3B+l+25^ where S denotes the spin, B the baryon number and 
L the lepton number of the particles. Hence Rp is equal to 1 for particles and equal to —1 for 
sparticles. We consider here the SUSY phenomenology at HERA in the presence of ^p Yukawa 
couplings but maintain otherwise the minimal field content of the Minimal Super symmetric 
Standard Model (MSSM) Of special interest for HERA are the Yukawa couplings between 
a squark and a lepton-quark pair . These are described in the superpotential by the terms 
X'^ji^LiQjDk, with i, j, k being generation indices]]. The corresponding part of the Lagrangian, 
expanded in fields, reads as: 



^LiQjDk - Kjk 



+ C.C. (1) 



where the superscripts c denote the charge conjugate spinors and the * the complex conjugate 
of scalar fields. Hence the couplings A'^^ allow for resonant production of squarks at HERA 
through eg fusion. For the nine possible couplings, the corresponding single production 
processes are given in table [l]. With an e+ beam, HERA is most sensitive to couplings X'lji, 
where mainly tt]^ squarks are being produced with a cross-section approximately scaling as 

'in the usual superfield notation, Li, Qj and Dk contain respectively the left-handed leptons, the left-handed 
quarks and the right-handed down quark, together with their SUSY partners l]^, q'j^ and d^. 



Table 1 : The two resonant squark pro- 
duction processes at HERA (e~^ beam) 
allowed by each i?-parity violating cou- 
pling a;^.^. 



X'lji ■ d{x) where d{x) is the probability to find a d quark in the proton with a momentum 

fraction x = M? /s and Mg denotes the squark mass. The production of the antisquark is 
also possible albeit with a much lower cross-section since a antiquark must participate in the 
fusion. 

The search presented here is performed under the simplifying assumption that one of the 
X[ji^ dominates. The squarks decay either via their Yukawa coupling into SM fermions (J^p), or 
via their usual gauge couplings (gauge decay) into a gluino g (the SUSY partner of the gluon), a 
neutralino Xa (« = 1, 4) or a chargino (/5 = 1, 2). The mass eigenstates Xa mixed states 
of the photino, the zino and the neutral higgsinos, which are the SUSY partners of the photon, 
of the Z and of the two neutral Higgs fields respectively. The charginos are mixed states 
of the charged higgsinos and of the winos, SUSY partners of the W^. Neutralinos, charginos 
and gluinos are unstable. This holds in SUSY also for the Lightest Supersymmetric Particle 
(LSP), assumed here to be a x (x° or x"*^) or a g, which decays into a quark, an antiquark and a 
lepton via a virtual squark or slepton undergoing a^p decay through the A' coupling. This is 
in contrast to Rp conserving SUSY models and has important phenomenological consequences. 

In cases where both production and decay occur through a A'^^^ coupling (e.g. Fig. [l|a and 
c for 7^ 0), the squarks have the same signature as scalar leptoquarks (LQ) [^. As can 

be seen from equation (|T]), the rf^ can decay either into e+ + or + d^, while the S-^ only 
decays into e+ + d''. The final state signatures consist of a lepton and a jet and are, event-by- 
event, indistinguishable from SM neutral current (NC) and charged current (CC) Deep-Inelastic 
Scattering (DIS). 

When the tt-^ (rf^) undergoes a gauge decay into a x°5 a x"*" or a ^ (a x° or a g) as shown 
in Fig. |T]b and d, the final state will depend on their subsequent decays. Neutralinos Xa ^^'^h 
a > 1 as well as charginos (gluinos) usually undergo gauge decays into a lighter x and two SM 
fermions (two quarks), through a real or virtual boson or sfermion (squark). The decay chain 
ends with the decay of one sparticle, usually that of the LSP. 

^.p decays of x's or gluinos are mainly relevant for the lightest states. Neutralinos can 
undergo the ^p decays x" ~^ e^qq' or x° i^qq, the former (latter) being more frequent if the 
X° is dominated by its photino (zino) component. Gluinos can undergo the same ]^p decays. 
When a x" or a ^ decays via ]^p into a charged lepton, both the "right" and the "wrong" charge 
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Figure 1 : Lowest order s-channel diagrams for I^p squark production at HERA followed by (a), (c) 
I^p decays and (b), (d) gauge decays of the squark. In (b) and (d), the emerging neutralino, chargino or 
gluino might subsequently undergo a Rp violating or Rp conserving decay of which examples are shown 
in the dashed boxes for (b) the Xi and (d) the Xi- 



lepton (with respect to the incident beam) are equally probable, this latter case leading to largely 
background free striking signatures for lepton number violation. In contrast, the only possible 
]^p decays for charginos are uu^d^ and ~^ e^dJ^dK 

The decay chains of tt-^ and c/^, analysed in this paper are classified into seven distinguish- 
able event topologies as described in table 0. This classification relies on the number of charged 
leptons and/or jets in the final state, and on the presence of missing energy. Channels labelled 
LQe and LQw are the "leptoquark-like" decay modes of the squark, proceeding directly via 
|?p, while the remaining channels cover the gauge decays of the squark and are characterised by 
multijet (MJ) final states. Channels labelled eMJ, e~MJ and uMJ involve one or two SUSY 
fermions (x or 9) denoted by X and Y in table ^ Channels e£MJ and ulMJ necessarily in- 
volve two SUSY fermions. Decay patterns involving more than two x or g are kinematically 
suppressed and are not searched for explicitly. The relative contributions of the channels consid- 
ered depend in particular on the value of the Yukawa coupling A' and on the gauginoQ-higgsino 
mixture of neutralinos and charginos. They will be shown as functions of the squark mass in 
section for some example cases. 



^The gauginos are the SUSY partners of the gauge bosons. 



Additional event topologies not listed in table |^ could in principle arise in the case where 
the Xi has such a small decay width (e.g. when it has large higgsino components) that it decays 
far away from the interaction point or leads to final states with displaced vertices . However 
the region of MSSM parameter space which would allow a to escape detection for a finite 
value of the coupling is now very severely constrained by the searches for charginos carried 
out at LEP p. The lifetimes of the sparticles are neglected in this analysis. 



3 The HI Detector 

A detailed description of the HI detector can be found in [^. Here we describe only the com- 
ponents relevant for the present analysis in which the final state of the events involves either a 
positron]] with high transverse energy or a large amount of hadronic transverse energy flow. 

The positron energy and angle are measured in a liquid argon (LAr) sampling calorime- 
ter [[1^] covering the polar angular range 4° < < 154° and all azimuthal angles. Polar angles 
are defined by taking the origin of the coordinate system to be at the nominal interaction point 
and the 2;-axis in the direction of the proton beam. The granularity of the LAr calorimeter is 
optimised to provide fine and approximately uniform segmentation in laboratory pseudorapid- 
ity 1] and azimuthal angle (p. The calorimeter consists of a lead/argon electromagnetic section 
followed by a stainless steel/argon hadronic section. Test beam measurements [ JTTl ] of the LAr 
calorimeter modules have shown an energy resolution of a{E)/E ~ 12%/^ E/ GeV © 1% 
for electrons and a{E)/E ~ 50%/^ E/ GeV © 2% for pions. The angular resolution on the 
positron measured from the electromagnetic shower in the calorimeter varies from ~ 2 mrad 
below 30° to < 5 mrad at larger angles. For the acquisition of events we rely on the LAr trig- 
ger system [ ]T^ ] whose efficiency is close to 100% for the transverse energies (Ex) considered 
here. A lead/scintillating-fibre backward calorimeter [|12|| extends the coverageQat larger angles 
(153° <e< 178°). 

The tracking system which is surrounded by the calorimeters is used in particular to deter- 
mine the position of the interaction vertex. The main components of this system are central drift 
and proportional chambers (25° < 9 < 155°), a forward track detector (1° < 6 < 25°) and a 
backward drift chamber. The tracking chambers and calorimeters are surrounded by a super- 
conducting solenoid providing a uniform field of 1. 15 T parallel to the z axis within the detector 
volume. The instrumented iron return yoke surrounding this solenoid is used to measure leak- 
age of hadronic showers and to recognise muons. In the very forward region {9 < 15°) muons 
can also be detected in three double layers of drift chambers, forming the Forward Muon De- 
tector. The luminosity is determined from the rate of Bethe-Heitler ep — * ep7 bremsstrahlung 
events measured in a luminosity monitor. 



4 Monte Carlo Event Generation 

For each possible SM background source, complete Monte Carlo simulations of the HI detector 
response are performed. Most of them correspond to a luminosity of more than 10 times that of 
the data. 

^Unless otherwise stated, the analysis does not distinguish explicitly between and e". 
^The detectors in the backward region were upgraded in 1995 by the replacement of the lead/scintillator tile 
calorimeter [ p3[ ] and a proportional chamber. 
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^ uqq 
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^ fi^fi^ Y 

A' 

"—>■ uqq 


£ (e or /i) 
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Table 2: Squark decay channels in SUSY classified by distinguishable event topologies. X and Y 
denote a neutralino, a chargino or a gluino. Quarks are generically denoted by q, except for the LQu 
channel which involves specific (s)quark flavours. The final states corresponding to ^ = r for the e£MJ 
and utMJ channels are not explicitly looked for in this analysis. 



For the simulation of the NC and CC DIS backgrounds, the DJANGO [|T^ event generator is 
used, which includes first order QED radiative corrections. QCD radiation is treated following 
the approach of the Colour Dipole Model [|T5p and is implemented using ARIADNE []T^. The 



hadronic final state is generated using the string fragmentation model [|T7]]. The parton densities 
in the proton used to estimate DIS expectations are taken from the MRST [ [TSl ] parametrisation. 



For direct and resolved photoproduction (7^) of light and heavy flavours, the PYTHIA event 
generator [[T^ is used which relies on first order QCD matrix elements and uses leading-log 



parton showers and string fragmentation [|T7|]. The GRV (GRV-G) parton densities in the 
proton (photon) are used. 

The simulation of the leptoquark-like signatures (LQe and LQu) relies on the event gen- 
erator LEGO pT| ] which is described in more detail in . For squarks undergoing gauge 



decays, we use the SUSYGEN [ ]23| ] event generator, recently extended to allow the gener- 
ation of SUSY events in ep collisions. Any gauge decay of the squark can be generated, and the 
cascade decays of the subsequent x's or g are performed according to the corresponding matrix 
elements. 

In both LEGO and SUSYGEN, initial and final state parton showers are simulated following 
the DGLAP [E5]] evolution equations, and string fragmentation [[TP|, [T7[] is used for the non- 



perturbative part of the hadronisation. In addition initial state bremsstrahlung in the coUinear 
approximation is simulated in the LEGO generator. The parton densities used [ p^ are evaluated 
at the scale of the squark mass. This scale is also chosen for the maximum virtuality of parton 
showers initiated by a quark coming from the squark decay. Moreover, in the SUSYGEN 
generator, the parton showers modelling QCD radiation off quarks emerging from a x or ^ 
decay are started at a scale given by the mass of this sparticle. 

To allow a model independent interpretation of the results, the signal topologies given in 
table 1^ were simulated for a wide range of masses of the SUSY particles. The events were 
passed through a complete simulation of the HI detector. The squark mass was varied from 
75 GeV to 275 GeV in steps of typically 25 GeV. Gauge decays of squarks involving one or 
two SUSY fermions (x or g) were simulated separately. For gauge decays of squarks into a a 
X"*" or a ^ which directly decays via_^p (i.e. processes corresponding to the first line of the eMJ, 
e~MJ and uMJ rows in table 0) the process q qxi was simulated for Xi niasses ranging 
between 40 GeV and 160 GeV. In order to study gauge decays involving two x or g, the process 
Q ~^ QXi ~* QXiff' was simulated for xt niasses ranging between 90 GeV and ~ M^, and 
for Xi niasses between half of the xt mass and ~ Masses of the x's were varied in steps 
of about 20 GeV. These simulations allowed the determination of signal selection efficiencies 
as a function of the masses of the squark and of the involved x 9 for essentially all allowed 
scenarios, since the grid size chosen for the simulated scenarios was small enough for a linear 
interpolation between them. 



5 Event Selection and Comparison with Standard Model 
Expectation 

The data reduction starts by the rejection of non-ep background, which is common to all chan- 
nels presented below. It is required that the events are accepted by a set of beam halo and cosmic 
muon filters p6|], that they satisfy constraints on their timing relative to the nominal time of the 
beam bunch crossings, and that a primary interaction vertex is reconstructed. 

Events containing lepton(s), hard jets, or a large amount of missing transverse energy are 
then selected using the following identification criteria: 



a positron (or electron) is identified by a shower shape analysis of clustered energy de- 
posits in the LAr calorimeter; the positron energy cluster should contain more than 98% 
of the LAr energy found within a pseudorapidity-azimuthal cone centered around the 
positron direction and of opening a/ (Ar/)^ + (A</))2 = 0.25, where rj = — lntan|; at 
least one charged track is required within this isolation cone; 

a muon candidate is identified as a track measured in the central and/or forward tracking 
system, which has to match geometrically an energy deposit in the LAr calorimeter com- 
patible with a minimum ionising particle, and/or a track in the instrumented iron and/or a 
track in the forward muon detector; 

hadronic jets are reconstructed from energy deposits in the LAr calorimeter using a cone 
algorithm in the laboratory frame with a radius ^/Arf~+~A(j? = 1; the fraction of the jet 
energy deposited in the hadronic part of the calorimeter must be at least 5%; 

the missing transverse momentum Pr.miss is obtained as 



PT,7niss = Y Ei sin Oi COS 0i j + Ei sin 6*^ sin cj)^ (2) 
where the summation runs over all energy deposits i in the calorimeters. 
In addition, the selection makes use of the following kinematic variables: 

• the momentum balance with respect to the direction of the incident positron, obtained as: 

Y,{E-P,) = Y, (1 - COS 9,) (3) 

where the summation runs over all energy deposits i in the calorimeters. J2 (E — Pz) 
should peak at twice the energy E'^ of the incident positron for events where only particles 
escaping in the proton direction remain undetected; 

• the Lorentz invariants y, and x characterising the kinematics of a DIS reaction, as well 
as the energy M in the centre of mass of the hard subprocess, are determined using the 
measurement of the polar angle 9e, the energy Ee and the transverse energy Et^c of the 
highest Et positron: 

_ _ i^e(l-COSge) ^2 _ _ Ol M - ^ 



the variables y,Q'^, x and M calculated using the Jacquet-Blondel ansatz [ETl] 



where Pt^h and J2iE — Pz)h calculated as in equations (0) and (Q), but restricting the 
summations to all measured hadronic final state energy deposits. 

The search for squarks decaying via^p couplings into channels LQe and LQu is identical 
to the search for first generation leptoquarks presented in Gauge decay channels are grouped 
into two classes, e + jets + X and u + jets + X. Preselection criteria are designed for these 
two classes of events, on top of which dedicated cuts are applied for the gauge decay channels 
listed in table ^. For all considered channels, the selection criteria are given in table ^ together 
with the resulting signal efficiencies and the numbers of observed and expected events. 



5.1 Analysis of Squark i?-Parity Violating Decays 



Channel LQe: Squarks decaying into the channel LQe have the same signature as scalar 
leptoquarks and are characterised by high NC DIS-like topologies. Such a process should 
manifest itself as a resonance in the measured Afg distribution, with a resolution of 3 to 6 GeV 
depending on the squark mass. The selection criteria are those described in [|l]]. The observed 
and expected mass spectra are shown in Fig. ^ to be in good agreement, with nevertheless a 
slight excess around 200 GeV already reported in [|l], The sources of systematic errors are 
described in section The (arbitrarily normalised) mass distribution expected from signal 
events coming from a 200 GeV squark decaying into the channel LQe is also shown. The peak 
value is slightly below the nominal squark mass due to final state QCD radiation Similar 
searches have been performed by the ZEUS experiment [p9|] . 
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Figure 2: Mass spectra for (a) NC DIS-like and (b) CC DIS-like final states for data (symbols) and DIS 
expectation (solid histograms). In (a) the comparison is shown after a Mg dependent cut on ye designed 
to maximise the significance of a squark signal [[T]]. The grey boxes indicate the it la band of systematic 
errors of the DIS expectations. The dashed histograms show the mass distributions for simulated events 
coming from a 200 GeV squark decaying into the channels (a) LQe and (b) LQu, with an arbitrary 
normalisation. 

Channel LQu: Squarks undergoing a LQv decay lead to CC DIS-like events with high 
missing transverse momentum showing a clustering in the Mh distribution with a resolution of 

about 10% of the squark mass. Only squarks, produced via a fusion between the incident 
e+ and a quark can undergo such a decay. The search for squarks decaying into the channel 
LQv is described in [jl]]. The observed and expected mass spectra are shown to be in good 
agreement in Fig. 0b. The (arbitrarily normalised) mass distribution expected from signal events 
coming from a 200 GeV squark decaying into the channel LQu is also shown. Similar searches 
have been performed by the ZEUS experiment p(]|]. 



5.2 Analysis of Squark Gauge Decays leading to e + jets + X Topologies 



When the squark undergoes a gauge decay leading to a positron, the final states can be clas- 
sified into several topologies, namely eMJ, e~MJ, eeMJ, ejxMJ and ueMJ. The "e- 
preselection" requirements which are common to all these e + multijet channels are the follow- 
ing: 
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Figure 3: Correlation between ye 
and the angle of the most backward 
jet, Obackward, for (a) SUSY Monte 
Carlo events where the squark un- 
dergoes a gauge decay leading to e -i- 
multijets +X final states and (b) NC 
DIS Monte Carlo events, when pre- 
selection cuts (1) to (3) are applied. 
Cut (4) only retains events above the 
diagonal line. In (a) events were 
generated for the range of masses 
considered in this analysis. 



1. at least one positron candidate in the angular range 5° < < 110° with E'^ e > 5 GeV; 

2. at least two jets in the angular range 7° < 9 < 145°; the highest Et jet must satisfy 
9jeti > 10° scad ETjeti > 15 GeV; thesecondhighest£'Tjetmusthave£JT,jet2 > 10 GeV; 

3. Ve > 0.4; 

4. of the two highest Et jets, the one with the larger polar angle, Obackward', must satisfy: 

Ve - 0.4 > {9backv.ard ' 25°)/60°; 

5. the minimum of the polar angles of the highest Et positron and of the two highest Et 
jets must satisfy: 

Min(^e,^jetl,^jet2) <45° . 

In gauge decays of a squark, a positron can emerge from the decay of a x or ^ appearing 
in the decay chain. It takes away a (possibly small) fraction of the momentum of this fermion, 
which motivates the cut (3). Moreover, it is strongly boosted in the direction of the incident pro- 
ton, such that the 9^ < 110° requirement discriminates the signal from the NC DIS background. 
Cut (4) exploits the fact that for high ye NC DIS events satisfying the above set of cuts, one 



hard jet is usually scattered in the backward region of the calorimeter. In contrast, jets coming 
from a squark gauge decay will be boosted in the forward direction. The Obackward distribution 
for SUSY events depends on the masses of the sparticles involved and cut (4) was designed to 
always retain a large fraction of the signal events. The effect of cut (4) is illustrated in Fig. |3|. 
Cut (5) requires that one of the squark decay products should be emitted in the forward direc- 
tion and allows an additional reduction of the SM background by ~ 40%, with a negligible 
efficiency loss on the signal. 
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Figure 4: For the e + multijets +X preselection, distributions of (a) the transverse energy Er^e of 
the highest Et positron; (b) Ue', (c) the transverse energy of the most forward jet; (d) the transverse 
energy of the most backward jet; (e) the polar angle of the most backward jet; (f) the missing transverse 
momentum. Superimposed on the data points (symbols) aie histograms of the SM expectation (DIS and 
7p). The grey band indicates the uncertainty on the SM prediction. The contribution from processes 
alone is shown as dotted histograms. 

Applying the above selection criteria, 214 events are accepted, which is in good agreement 
with the SM prediction of 210 ± 34, including 47 events from photoproduction where a jet 
has been misidentified as an electron. Fig. ^ shows the observed distributions of the transverse 
energy of the highest positron, y^,, the transverse energy of the most forward and most 
backward jet, the polar angle of the most backward jet and the missing transverse momentum. 
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Figure 5: Mass spectrum for e + multijet topology eMJ for the data (symbols) and the SM expec- 
tation (solid histogram). The grey band indicates the uncertainty on the SM prediction. The dashed 
histogram shows the expected mass distribution for events coming from a 200 GeV squark decaying into 
the channel eMJ, with an arbitrary normalisation. 

All distributions are seen to agree well with the SM expectation within the systematic errors. 

For channels leading to e + jets + X final states additional cuts, listed in table ^ are 
applied on top of the preselection requirements (1) to (5). In each case, good agreement is ob- 
served between the data and the SM expectation largely dominated by the NC DIS contribution. 
Additional information for the different channels is given below. 

Channel eMJ: A mass Mj„^ is calculated as: 



where the sum runs over all energy deposits in the calorimeters for 9 > 10°, thereby excluding 
the proton remnant. For squarks decaying into the eMJ channel, Mj„y provides an estimate of 
the q mass. This reconstruction method yields a typical resolution of 7 to 10 GeV depending 
on the squark mass. The Miny spectrum of the selected events is shown in Fig. |^ to be well 
described by the SM prediction. Also shown is the (arbitrarily normalised) mass distribution 
expected from signal events coming from a 200 GeV squark decaying into the channel eMJ. 
No charge determination of the lepton is performed here. 

Channel e~MJ: We consider the track in the e isolation cone which has the highest momen- 
tum projected on the axis defined by the event vertex and the centre of gravity of the calorimetric 
energy deposits associated with the electron. This track is required to have a reliably measured 
negative charge. The efficiency of the track quality requirements is ~ 80%, derived from data 
(candidates for channel eMJ ) and well reproduced by the NC DIS simulation. 

Channel e£MJ: Di-lepton final states are searched for provided that the lepton accompany- 
ing the e belongs to the first or second generation. For the channel efiMJ, the properties of 
the muons observed in the preselected events over the full range in transverse momentum were 
found to be well described by the simulation, as exemplified in Fig. ^. 
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Figure 6: Distribution of the transverse 
momentum Pt^^ of the muon for events sat- 
isfying the preselection criteria and where a 
muon has been identified in the angular range 
10° < 6* < 110°, for data (points) and SM 
(histogram). 



Channel veMJ: The common preselection criteria are complemented by a PT,miss require- 
ment and by the cut ?/e ■ (?/e — > 0.05. A cut on the product of ?/e with the difference {ye — Vh) 
exploits the fact that, for events coming from a squark decaying into the channel veMJ, the 
escaping neutrino carries a non-negligible part ofJ^i^^Pz) and hence the variable yh is sub- 
stantially smaller than y^, while ye ~ yn is expected for NC DIS events. Fig. ^ shows the 
distribution of ye ■ {ye — yh) for the 214 events accepted by the preselection and for the SM 
expectation. 
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Figure 7: Distribution of the variable ye • (ye — Vh) for the 214 events satisfying the preselection criteria 
(symbols) and for the expectation from NC DIS and 7p processes (open histogram). The grey error 
band on the open histogram indicates the uncertainty on the SM prediction. The arbitrarily normalised 
hatched histogram shows how this variable is distributed for events coming from a squark decaying into 
the channel ueMJ, generated for the range of masses considered in this analysis. 



5.3 Analysis of Squark Gauge Decays leading to v + jets + X Topologies 

Two channels are considered to cover cases where squarks undergo a gauge decay leading to a 
neutrino (and no positron) in the final state. These final states, uMJ and ujxMJ, are selected 
by the following 'V-preselection" requirements: 

1. a missing transverse momentum PT,miss > 25 GeV; 

2. at least two jets in the angular range 7° < 9 < 145° and with Et > 10 GeV, with the 
highest £'t jet satisfying Ojeti > 10° and Erjeti > 15 GeV. 
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Figure 8: For the u + multijets +X preselection, observed (symbols) and expected (open histograms) 
distributions of (a), (b) the transverse energies of the two highest Et jets and (c) the polar angle of the 
highest Et jet. The hatched histograms show distributions for events from a squark gauge decay into 
u + jets, generated for the range of masses considered in this analysis, (d): The mass Mj-ec, corresponding 
to the energy in the centre of mass of the hard subprocess assuming that only one neutrino escapes 
detection; the dashed histogram shows the Mrec distribution for signal events coming from the decay of 
a 200 GeV squark into the channel uMJ. All histograms showing the SUSY expectation are arbitrarily 
normalised. 



We observe 44 events satisfying these preselection criteria, in good agreement with the SM 
prediction of 46.5 ± 6.9, mainly coming from CC DIS processes. 

Figs. ||a-c show the distributions of the transverse energies of the two highest Et jets and of 
the polar angle of the highest Et jet. The data are well described by the SM expectation. Dis- 
tributions for squarks decaying into u + jets are also shown. Assuming that the missing energy 
is carried by one neutrino only, its kinematics is reconstructed exploiting energy-momentum 
conservation. The four- vector of this u is then added to that of the hadronic final state to recon- 
struct the invariant mass Mrec of the incoming electron and quark. For squarks decaying into 
the channel uMJ , Mrec provides an estimate of the squark mass. The observed and expected 
distributions for Mrec are in good agreement as shown in Fig. |8|d. The dashed, arbitrarily nor- 
malised, histogram in Fig. ||d shows the resulting mass spectrum for a 200 GeV squark decay 
into the uMJ channel. The observed resolution of ~ 15 GeV is typical for the range of squark 
masses probed in this analysis. 

Final cuts and results for the channels wMJ and w^MJ are given in table ^. The number 
of observed candidates is in good agreement with the SM expectation, which is dominated by 
the contribution of CC DIS processes. 



5.4 Systematic Errors 

In each channel, the error on the expectation from Standard Model processes has been calculated 
by taking into account the systematic errors described below. The experimental error sources 
considered are: 



an uncertainty of ±1.5% on the integrated luminosity; 

an uncertainty on the absolute calibration of the calorimeters for electromagnetic ener- 
gies, ranging between ±0.7% in the central part and ±3% in the forward region of the 
LAr calorimeter; this constitutes the main error source for the estimation of NC DIS 
background to the channel LQe; 

an uncertainty of ±4% on the absolute hadronic energy scale. For inclusive NC DIS final 



states, the over-constrained kinematics allows a reduction of this uncertainty to ±2% [[2^], 
which applies to channels LQe and LQw. This is the main error source for all channels 
except LQe. 

The following theoretical uncertainties on SM cross-sections are considered. 

• For NC DIS-like final states, an uncertainty of ±5% is attributed to the proton struc- 



ture [pSQ, which is partly due to the experimental errors on the input data entering the 
QCD fits, and partly linked to the assumptions for the shapes of the parton distributions 
at the scale where the perturbative QCD evolution is started. For CC DIS-like topologies, 
which are mainly induced by d quarks whose density in the proton is less constrained, this 
uncertainty increases linearly with up to ~ 20% at the highest Q"^ considered here [jT]]. 
In adition, the error on the strong coupling constant as leads to an uncertainty of ±4% 
on the proton structure. This was inferred [EHp by comparing the CTEQ4 (Al) to (A5) 



parametrisations [ |3T1 ] with as{Mz) ranging between 0.110 to 0.122. 



Channel 


Selection Cuts 


Efficiency 


Nobs 


N 


DIS-like channels: Q'^ > 2500 GeV^ y < 0.9 


LQe 


ET,e > 15 GeV 

PT,miss/'^jET,e < W GeV 

40 < ^ - P,) < 70 GeV 
optimised lower y-cut 


40-70 % 


310 


301 ± 23 


LQu 


PT,m^ss > 30 GcV 

no electron E^^e > 5 GeV 


30-80 % 


213 


199 ± 12 


channels with: e + multijets + X 

e-preselection: E^^e > 5 GeV; > 2 jets: E'rjeti,2 > 15, 10 GeV; high y^,; angular cuts 


eMJ 


PT,miss < 20 GeV 
40 < E - Pz) < 70 GeV 


35-50 % 


159 


151 ± 17 


e-MJ 


eMJ criteria 
+ "wrong" charge of e 


^ 30% 





1.3 ±0.5 


eeMJ 


second e with: 

ET,e2 > 5 GeV 
5° < 9^2 < 110° 


^ 30% 





0.7 ± 0.4 


e^MJ 


Pt,^ > 5 GeV 
10° <9^< 110° 


35-50% 


2 


4.2 ± 1.2 


veMJ 


PT,miss > 15 GeV 
yeiye - yh) > 0.05 


^ 30% 


1 


3.2 ± 1.2 


channels with: u + multijets + X 

z/-preselection: PT,m,iss > 25 GeV; > 2 jets: ETjeti,2 > 15, 10 GeV 


vMJ 


ET,jet2 > 15 GeV 
EiE- P,)^ < 55 GeV 


20-60 % 


21 


23 ±4 




Pt,^ > 5 GeV 
10° <9^< 110° 


^40% 





0.5 ± 0.2 



Table 3: Selection criteria, typical efficiencies, total number of observed events Nobs and the cor- 
responding SM expectation N^xp with its uncertainty, for each squark decay channel analysed. The 



"e-preselection" criteria are detailed in section 5.2. 



• Higher order QED corrections imply a ±2% uncertainty in the y range considered here [^Bp, 
estimated using the HECTOR [ p2| ] program. 



• An uncertainty of ±10% on the predicted cross-section for multijet final states is esti- 
mated by comparing leading order (LO) Monte Carlo simulations where higher order 
QCD radiation is modelled by either the dipole model or DGLAP parton showers, and 
next-to-leading order (NLO) calculations. 

For each error source, the analysis has been repeated shifting the central value by ±1 standard 
deviation to estimate their individual contribution. The overall systematic error on SM expec- 
tations is then determined as the quadratic sum of these individual errors and the statistical 
uncertainty on the Monte Carlo simulation. 



6 Constraints on SUSY Models 

No significant deviation from SM expectations has been found in the analysis of various I^p and 
gauge decay channels. These channels are combined to set constraints on I^p SUSY models. 

As mentioned in section ^, when HERA operates with incident positrons, the best probed 
couplings are A'^^^, on which we mainly concentrate here. Such a coupling allows for the 

production of u'^ or rf^ squarks, the latter with a much reduced rate due to the smaller parton 
density. Thus only u'j^ production is considered. Since this squark does not decay into v + q the 
channel LQu is not taken into account in the derivation of limits. The sr production is also 
(conservatively) neglected when setting limits on A'^g- 

In this section, mass dependent upper limits on the Yukawa couplings are first derived in a 
"phenomenological" version of the MSSM where the masses of the sfermions are not related to 
the SUSY soft-breaking mass terms of the gauginos. Scans are then performed in the framework 
of this "phenomenological" MSSM as well as in the constrained MSSM, and bounds on the 
Yukawa couplings are set in these models. Results are finally interpreted in the framework of 
the minimal Supergravity model. 

6.1 Derivation of Limits 

Mass dependent upper limits on the production cross-section a{e^p u-'j^) are obtained as- 
suming Poisson distributions for the SM background expectations as well as for the signal. A 
standard Bayesian prescription with a flat prior probability density for the signal cross-section is 
used. Both systematic and statistical errors have been folded in channel by channel as described 
in [0]. Each channel contributes in the derivation of the limits via its branching ratio, the num- 
ber of observed and expected events satisfying the relevant selection cuts and the corresponding 
selection efficiency. For the channels LQe, eMJ and vMJ, the numbers of observed and 
expected events are integrated within a mass bin which slides over the accessible mass range. 
The width of the mass bin is adapted to the expected mass resolution in each channel, such that 
this bin contains approximately 68% of the signal at a given squark mass. For the channels 
e~MJ, eiMJ and ulMJ, where both the SM expectation and the observation are small, no 
mass restriction is imposed. 



Events which fulfil the selection requirements of more than one channel are only counted 
as squark candidates in the channel with the highest sensitivity. This prescription is illustrated 
in table 0. Note that the potential overlap between channels LQe and eMJ is already con- 
siderably reduced by the mass bin requirements, since the corresponding reconstructed squark 
masses, Me and Mj„„ respectively, differ by typically more than 20 GeV. The relative efficiency 
loss on the SUSY signal induced by this slight tightening of the cuts depends on the masses of 
the sparticles involved, and is found to vary between ~ 0.5% and ~ 3%. 



candidates fulfilling the selection criteria for channels: 


are not considered in channel: 


eMJ e-MJ eeMJ efiMJ veMJ 


LQe 


e'MJ eeMJ efiMJ ueMJ 


eMJ 


eeMJ efiMJ 


veMJ 


ly^MJ 


vMJ 



Table 4: Prescription adopted to ensure no overlap between the considered channels. 



The masses of the neutralinos, charginos and gluinos, as well as the couplings between 
any two SUSY particles and a standard fermion/boson, are determined by the usual MSSM 
parameters: the "mass" term jj, which mixes the Higgs superfields, the SUSY soft-breaking 
mass parameters Mi, M2 and M3 for f/(l), SU{2) and SU{?)) gauginos, and the ratio tan/5 
of the vacuum expectation values of the two neutral scalar Higgs fields. These parameters are 
defined at the electroweak (EW) scale. We assume that the gaugino mass terms unify at a Grand 
Unification (GUT) scale to a common value mi/2 leading to usual relations between Mi, 
M2 and M3, and approximate the gluino mass by the value of A-'L^ at the EW scale. The masses 
and decay widths of all involved sparticles have been obtained using the SUSYGEN package. 

For a fixed set of MSSM parameters and sfermion masses, the branching ratios for the dif- 
ferent channels only depend on the Yukawa coupling. This also holds for the upper limit on 
the signal cross-section derived from the combination of the analysed channels. At a given 
squark mass, Yukawa couplings for which is smaller than the signal cross-section are ex- 
cluded, where the signal cross-section is obtained by multiplying the leading-order production 



cross-section by i^-factors [ pBQ accounting for NLO QCD effects. These can enhance the signal 
rate by C(10%). 

Decay chains involving more than two SUSY fermions (x or g) can contribute in principle 
to the gauge channels analysed. In these cases parameterising the signal efficiencies is not 
straightforward. Hence, only cascades involving two SUSY fermions are taken into account in 
the calculation of the visible branching ratios for gauge decay channels. This determination of 
the branching ratios is conservative. It has been checked that the visible branching is generally 
close to 100%. Decays of x's into a Higgs boson are included in the calculation of the visible 
branching ratios when the Higgs decays into hadrons. The contribution of these decays is 
however very small. Hence the limits do not depend on the mass niA of the pseudoscalar Higgs, 
set here to 300 GeV in the models where is not related to the other parameters. 

The case of a non-vanishing coupling allowing for the resonant production of a stop 
squark (SUSY partner of the top quark) has to be treated separately. Firstly, the large top mass 
can not be neglected in the calculation of the branching ratios for the decays t or t ^ tg, 

which may eventually be kinematically forbidden. Secondly, the top quark decays via t — > hW . 
Most of the stop decays are in fact covered by our analysis, but the efficiencies for the considered 



channels, which are valid for decay patterns as shown in table ^ can not be used in that case. 
Conservatively, diagrams which lead to a top in the final state are thus not taken into account 
in the calculation of the visible branching ratios. For example, the decays Xi ~^ e~i:J will 
not be included in the branching ratios for the eMJ and e~MJ channels even when these 
are kinematically allowed. As a result, only the neutrino decays of the x'i (x'i ~^ ydd) will 
contribute in the derivation of limits on A'l 01 . 



6.2 Limits on A^^^ and A^^2 "phenomenological" MSSM 

We consider here a version of the MSSM where the parameters /i, M2 and tan (3 are only used to 
determine the masses and couplings of the x's, while the sfermion masses are free parameters. 
We neglect any possible mixing between sfermions and assume that all squarks are degenerate in 
mass. This assumption only enters in the calculation of the branching ratios of the x's and of the 
gluino, since we are mainly probing the tt^ squark. Sleptons are also assumed to be degenerate, 
and their mass Mj is set either to the common squark mass, or to a fixed value (90 GeV) close 
to the lowest mass bound from sfermion searches at LEP. We first derive constraints on the 
couplings y^'iji, where a squark could be produced via an e^d fusion, and consider in a second 
step squark production via e+s fusion through a coupling. 

Example upper limits obtained at 95% confidence level (CL) on A'^^^ (j = 1, 2) as a function 
of the u'l^ mass are shown in Fig. |^a and Fig. [T^a, under the assumption = 90 GeV or 



Mf = Mq. The MSSM parameters are chosen such that the lightest neutralino is dominated by 
its photino component in Fig. |^ and by its zino component in Fig. [1^. The gluino mass is large 
due to the large M2 (and hence M3) values and thus squark decays into g are kinematically 
forbidden. In the four scenarios considered, Yukawa couplings larger than ~ 0.01 (~ 0.5) are 
ruled out for squark masses of 75 GeV (275 GeV). 

The relative contributions of all channels are shown in Fig. |^,c and Fig. [T(]|b,c against the 
squark mass, at the current sensitivity limit on the Yukawa coupling. In each case the total 
branching ratio covered by the analysis is close to 100 %. At the largest squark masses, where 
a large Yukawa coupling is necessary to allow squark production, the relative contribution of 
the LQe channel becomes important. In the case illustrated in Fig. ^ (Fig. [1^ the dominant 
channels are those with an (z/) in the final state, resulting in particular from the main decay 
mode of the Xi- The relative contributions of the channels strongly depend on the slepton mass. 
In the case illustrated in Fig. where a small slepton mass = 90 GeV is used, the two- 
body decay of the ~ 95 GeV xt or X2 ii^to ^ lep ton- slepton pair is kinematically allowed. As 
a result, the contributions of the channels e£MJ and uiMJ are considerably enhanced. In 



the case shown in Fig. |T0|, the channel uMJ remains dominant even in the light sleptons case 
because the xt mass is here ~ 80 GeV. Only squark decays into the heavier x^ lead to the 
enhancement of the channel ulMJ shown in Fig. |TD|c. 

Despite the fact that the relative contributions of the various channels are strongly model 
dependent, the upper limits on the Yukawa coupling do not depend significantly on the scenario 
considered because the sensitivity of our analysis is similar in all gauge channels, and because 
the covered branching ratio is always close to 100%. 

In order to investigate more systematically how the sensitivity depends on the MSSM pa- 
rameters, a scan of the parameters M2 and fi is performed, for tan/3 = 2. The effect of 
varying the parameter tan (3 will be studied in the next section. The mass of the sleptons is 



set to 90 GeV, the parameters M2 and fi are varied in the range 70 GeV < M2 < 300 GeV and 
-300 GeV < /i < 300 GeV. Points which lead to a scalar LSP or to LSP masses below 30 GeV 
are not considered. This latter restriction, as well as the lower value for M2, are motivated by 
the exclusion domains resulting from x searches in SUSY at LEP. For each point in this 
(fi, M2) plane the upper bound X'^^ on the coupling A'^^^-^ is obtained. The results are shown in 
Fig.[TT]a for A'^^^^ (j = 1, 2) and in Fig.pTjb for X'l^i- The two full curves in Fig.|Tl] indicate the 
maximal and minimal values obtained for AJ^^ within the parameter space investigated. 

The spread between these extrema for X'^^ is small for squark masses above 150 GeV and 
decreases with increasing squark mass. Comparing Fig.pT|a and FigjTT|b, the constraints on 
X'l^i and on A'^^^^^ (7 = 1,2) are seen to be quite similar. Only for small squark masses is the 
sensitivity on A'^^g^^ reduced because of the small efficiency in the wMJ channel. 
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Figure 9: (a): Upper limits at 95% CL for the coupling X'^^^ (j = 1, 2) as a function of the squark 
mass, for a set of MSSM parameters leading to a x? of ~ 40 GeV dominated by its photino component. 
Regions above the curves are excluded. The limits are given for two hypotheses on the slepton mass, 
(b): The branching ratios of all channels versus the squark mass, when sleptons and squarks are assumed 
to be degenerate; (c): as (b) but assuming a slepton mass of 90 GeV. The thick curves in (b) and (c) 
indicate the total branching covered by the analysis. 



For a Yukawa coupling of electromagnetic strength Q^^-^/Att = Uem, i-e- ^'iji = 0.3) squark 
masses below ~ 260 GeV are excluded at 95% CL . This extends beyond the mass domain 
excluded from relevant searches for scalar leptoquarks performed by the DO CT] and CDF p5|] 



experiments, which rule out u'^ squark masses below 205 GeV if the branching ratio of the 
squark into eg is greater than ^ 50%, and below ~ 110 GeV when this branching is ~ 10%. 
Since in SUSY such a branching can be naturally small (as seen above) such leptoquark- 
like constraints are not very stringent. Direct squark searches at LEP exclude masses below 
~ 90 GeV. Bounds from searches for SUSY carried out at the Tevatron do not apply in the 
unconstrained model considered here. 

Our results are also compared in Fig. [njto indirect limits [^. The production of a u squark 
via a X'lii coupling is very severely constrained by the non-observation of neutrinoless double 
beta decay P7|]. The best indirect limit on the coupling A'^gi C-^'isi)' which could allow for the 
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Figure 10: (a): Upper limits at 95% CL for the coupling A'^^^ {j = 1, 2) as a function of the squark 
mass, for a set of MSSM parameters leading to a of ~ 40 GeV dominated by its zino component. 
Regions above the curves are excluded. The limits are given for two hypotheses on the slepton mass, 
(b): The branching ratios of all channels versus the squark mass, when sleptons and squarks are assumed 
to be degenerate; (c): as (b) but assuming a slepton mass of 90 GeV. The thick curves in (b) and (c) 
indicate the total branching covered by the analysis. 



production of squarks c (t), comes from atomic parity violation (APV) measurements [j^, P^]. 



For squark masses below ~ 240 GeV the HI direct limits significantly improve this indirect 
constraint on \[2i by a factor up to ~ 3. In the case of a non- vanishing coupling, our 
results are more stringent than the APV constraints for stop masses between ~ 100 GeV and 
~ 220 GeV 

The limits on the couplings A'^^^ can be translated into upper bounds on the couplings \[j2^ 
which could allow for the production of a resonant via a fusion between the positron and a 
strange quark coming from the proton. To obtain the limit on \[j2 the ratio of the s and d quark 
densities and its uncertainty were taken from a LO QCD fit similar to [0^], in which various 



neutrino DIS data providing constraints on s{x) were considered. Resulting upper bounds at 
Mg = 200 GeV are shown in table ^. The limits were conservatively derived taking into account 
the 2cr uncertainties of the parton densities as given by the fit. Using the central prediction for 
s(x) as given by the MRST parametrisation, limits on \[j2 do not change by more than 20%. 
Table |] also shows existing indirect bounds [ ]3^ ] for comparison and recalls the bounds obtained 



on The sensitivity of our analysis to the coupling is significantly better than that 

coming from the leptonic decay width of the Z boson. No attempt was made to derive limits on 
couplings X[j.^ due to the large uncertainties on the b quark density at such high x and Q^. 
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Figure 11: Upper limits at 95% CL on the coupling (a) A'^^^^ (j=l»2) and (b) as a function of 
the squark mass for tan /3 = 2, in the "phenomenological" MSSM. For each squark mass, a scan of the 
MSSM parameters M2 and fi has been performed and the largest (lowest) value for the coupling limit 
is shown by the upper (lower) full curve. The dotted curve in (a) indicates the indirect bound on A'^^^^ 
from neutrinoless double beta decay assuming a gluino mass of 1 TeV. The dashed-dotted curves show 
the indirect bounds from atomic parity violation. 
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Table 5: 95% CL upper bounds on the couplings A'^^^ and A'^j2 ^ squark mass of 200 GeV. Also 
shown are the indirect bounds obtained from neutrinoless double beta decay, atomic parity violation, 
charged current universality, the upper bound on the neutrino mass, and the leptonic decay width of the 
Z boson. 



6.3 Limits on A^^^ in the Constrained MSSM 

In this section we consider a "constrained" (supergravity inspired) version of the MSSM where 
the number of free parameters is reduced by assuming, in addition to the GUT relation men- 
tioned in section ^TTJbetween Mi, M2 and M3, a universal mass parameter mo for all sfermions 
at the GUT scale. The evolution of the sfermions masses towards low scales is given by the 
Renormalisation Group Equations (RGE) and depends on the gauge quantum numbers of the 
sfermions. As a result, the sfermions masses at the electroweak scale are related to each other 
and to the parameters determining the gaugino sector. The model is thus completely deter- 
mined by e.g. mo, M2, n and tan (3 {tua is set to 300 GeV and we assume no mixing between 
the sfermions at the electroweak scale). 

For a given value of the £t]^ squark mass, the requirement of sfermion unification at large 
scale imposes an upper bound on the parameter M2, which is obtained using approximate so- 
lutions for the RGEQ. The upper bound on M2 increases with the squark mass and is smaller 
for the stop than for other squarks. As before sets of parameters leading to a scalar LSP or to 
a LSP mass below 30 GeV are not considered. This lower bound on the LSP mass forbids too 
small values of M2 and hence imposes a lower bound on the tt]^ mass, which is more stringent 
in case of the stop. A scan of /i, tan (3 and M2 is performed within —300 GeV < /i < 300 GeV, 
2 < tan /3 < 40 and within the M2 range allowed by tan (3 and the u'^ mass. 

The curves in Fig.[T^ indicate the maximal and minimal values for X[^^ as a function of 
the u'j^ mass. The spread of the domain spanned by the limits X'^^ is quite small, i.e. the 
sensitivity of our analysis on A' does not depend strongly on the free parameters of the model, 
in particular on tan (3. The most stringent limits are usually obtained for intermediate tan (3 and 
are in general better than those derived previously in the "unconstrained" MSSM because in this 
range the sneutrinos can be very light, leading to an enhancement of the quasi background-free 
channels e£MJ and vlMJ via e.g. Xi l^^- 

For a Yukawa coupling of electromagnetic strength, squark masses up to values of 260 — 
270 GeV can be ruled out at 95% CL in the framework of the constrained MSSM. Moreover, 
for a coupling strength 100 times smaller than a^.^, the most conservative bound on the mass of 
the u'^ obtained from the present analysis still reaches 182 GeV. 

Searches for SUSY performed at LEP [|]] also set limits on the model considered here. 
At LEP the mass domain explored by direct searches for squarks with couplings is limited 
by the beam energy. However searches for neutralinos and charginos lead to a lower bound on 
M2 which, using the RGE's, can be translated into a lower bound of ~ 240 GeV on the squark 
mass, thus reducing the allowed mass domain probed in Fig.[T^ where only the less stringent 

^ The possible effect of the Yukawa couplings A'j^ ^, on the RGE has not been taken into account here. 




Figure 12: Upper limits at 95% CL for the coupling A'^^^^ as a function of the mass of the tt]^, in the 
constrained MSSM. For each squark mass, a scan of the MSSM parameters fi, M2 and tan /3 has been 
performed and the largest (smallest) value for the coupling limit is shown by the upper (lower) curve. 
The resulting band is contained within the thin curves for A'^^^^^ (j = 1, 2) and within the thick ones for 
X'l^i- The requirement on the LSP mass imposes the u-'^ to be heavier than ~ 160 GeV (~ 235 GeV) for 
j = l,2(j = 3). 

condition Mlsp > 30 GeV was imposed. The combined searches for x's and sleptons at LEP 
moreover increase the lower M2 bound for low values of the mo parameter. The resulting lower 
bound on the mass of first and second generation squarks is close to the current HERA centre 
of mass energy. 



6.4 Constraints in the Minimal Supergravity Model 

The model considered above can be further constrained by imposing a common SUSY soft- 
breaking mass term for all scalar fields, and by assuming that the breaking of the electroweak 
symmetry is driven by radiative corrections. These additional assumptions lead to the so-called 



minimal supergravity (mSUGRA) model [^. By requiring Radiative Electroweak Symmetry 
Breaking (REWSB), the modulus of fi is related to the other model parameters. The program 
SUSPECT 1.2 [ p4l | is used to obtain the REWSB solution for when the other parameters are 
fixed. 

Assuming a fixed value for the coupling A'^^^, our searches can be expressed in terms of 
constraints on the mSUGRA parameters, for example on (mo, mi/2) when tan /3, the common 
trilinear coupling at the GUT scale Aq, and the sign of fx are fixed. The parameter Aq enters 



only marginally in the interpretation of the results and Aq is set to zero. Values of the param- 
eters leading to a LSP lighter than 30 GeV have not been excluded here. However a vanishing 
efficiency has been assumed for squarks undergoing a gauge decay ending by a x or ^ lighter 
than 30 GeV, since the parametrisation of the efficiencies (see section ^ is not valid in this 
domain. 

For tan/5 = 2 and /x < 0, results obtained for a Yukawa coupling A^^^ = 0.3 (j = 1, 2) 



are shown in the (mo, mi/2) plane in Fig. |13|a. The parameter space where M-j < 260 GeV 




mo (GeV) mo(GeV) 



Figure 13: Domain of the plane (mo,mx/2) excluded by this analysis for fi < 0, Aq = and (a) 
tan/3 = 2 or (b) tan/3 = 6, for a^p coupling A'^^^ = 0.3 (j = 1,2) (light shaded areas). In (b) the 
exclusion domain obtained for X[j^ = 0.1 (j = 1, 2) is also shown as the dark grey area. The hatched 
domains correspond to values of the parameters where no REWSB is possible or where the LSP is a 
sfermion. The regions below the dashed curves are excluded by the DO experiment; in (a) also the L3 
bound is shown as the dashed-dotted curve; these limits do not depend on the value of the coupling. 
Two isolines for the mass of the it-'^ are also shown as dotted curves. 

is nearly fully excluded. At low mi/2 values where the lightest x's and the g are lighter than 
30 GeV, the sensitivity on the m-^ mass however decreases since the efficiency is conservatively 
set to zero for all channels but LQe, and reaches ~ 210 GeV. Fig. [T3|a also shows the domain 
excluded by the DO experiment from searches for SUSY where Rp is violated by a A'^^^ 
coupling, relying on di-electron events. HI and Tevatron results are quite similar at low mo. 
However, the mSUGRA parameter space is still more constrained by the combined searches 
for x's and sleptons carried out at LEP as shown in Fig. [T3|a. LEP and Tevatron bounds do not 
depend on the value of the Yukawa coupling. 

Results for tan/3 = 6 and A'^^-,^ = 0.3 or 0.1 (j = 1,2) are shown in Fig. [r3|b. The ex- 
cluded domains extend considerably beyond the region ruled out by the DO experiment. This 
is due to the fact that, for large values of tan /3, the lightest neutralino is dominated by its zino 
component, so that its decay into is suppressed. As a result the sensitivity of the di-electron 
DO analysis is decreased, while the dominant squark decay mode is still observable in the HI 



analysis via the uMJ and uiMJ channels. LEP limits in [|p have not been given for this 
value of tan P but the corresponding bounds on mi/2 are expected to be similar to those shown 
in Fig. |T3|a within ~ 10 GeV. 

We now consider a non-vanishing coupling X[^^ which could lead to the production of a 
stop. The weak interaction eigenstates and mix in this case through an angle 9t to form 
two mass eigenstates, labelled ii and t2 (^i being the lightest by convention): 

ii \ ^ [ cos9t sin 6*4 \ f h 
t2 J \ - sin 6t cos 6t J \tii 

The production cross-section of the ii (t2) scales as Aj^l^^ cos^ 6*4 (A^l^ sin^ 6'^) since only ii 
enters in the LiQ^Di operator. Hence, the lightest stop ii does not necessarily have the largest 
production cross-section. Thus both ii and ^2 are searched for in the analysis. 

For channels e~MJ, eiMJ and ulMJ where the signal is integrated over the whole mass 
range the fraction of the visible signal in a given channel, k, is '^i^i2l3k,i£k,i(^i/ (^tou where 
is the branching ratio for tj to decay into this channel k, ek,i the corresponding selection 
efficiency, cTj the production cross-section of U, and atot = cri + cr2 the total signal cross-section. 

For the channels LQe, eMJ and uMJ where the signal is integrated over a "sliding mass 
bin" only the contribution of the state ti for which the sensitivity is maximal (i.e. which max- 
imises cri(^;, Pk,i£k,i)) is taken into account in the above summation. The numbers of observed 
and expected events are then integrated in the mass bin corresponding to ti only. 



The results are shown in Fig. IT^ for X[^^ = 0.3, Aq = 0, fi < and tan /? = 2 or tan f3 = Q. 



The domain below the line mi/2 ^ 10 GeV is not considered since it corresponds to cases 
where the only allowed LSP decay into ubd is kinematically forbidden. For tan (3 = 2, the 
excluded domain is slightly larger than that ruled out previously for A'^^^^^ = 0.3 (j = 2), due 
to the mixing in the stop sector which leads to ti masses smaller than the masses of the other tt^ 
squarks. In particular, larger values of niQ can be probed. This remains the case for tan /3 = 6 as 
long as mi/2 is large enough to ensure that the mass of the lightest neutralino is above 30 GeV. 
When the x? becomes too light, the efficiencies for the channels involving ^xt ^ yH decay 
are set to zero, and the sensitivity to the signal is only provided by the LQe channel or by the 
decays i bxt followed by a^p decay of the chargino. As a result, only lighter stops can be 
probed, for which the visible cross-section is large enough. Note that for both values of tan /5, 
masses of ii up to 245 GeV can be excluded for \[^^ = 0.3. This is slightly smaller than the 
maximal sensitivity of ~ 260 GeV reached, for the same coupling value, on the tt]^ = 1)2) 
mass (Fig. [T3|), or on the mass in the constrained MSSM when the stop mixing is neglected 
(Fig. O). This is due to the cos^ 9t reduction of the ii cross-section. 



For intermediate values of mo, the L3 sensitivity is slightly better than the limit obtained 
from this analysis for tan/3 = 2. In the same part of the parameter space, the HI limit is 
comparable with or slightly extends beyond the expected LEP sensitivity for tan P = 6. 



7 Conclusions 

We have searched for squarks in e~^p collisions at HERA in i?-parity violating SUSY models. 
No evidence for the resonant production of such particles was found in the various channels 
considered which cover almost all decay modes. Mass dependent limits on i?-parity violating 




Figure 14: Domain of the plane {mQ,mi/2) excluded for a. coupling = 0.3, for fi < 0, 
^0 = and (a) tan /5 = 2 or (b) tan f3 = 6. The hatched domains on the left correspond to values of 
the parameters where no REWSB is possible or where the LSP is a sfermion. For too small (the 
domain below the line mi^2 ~ 10 GeV), the LSP becomes lighter than the b quark and thus is stable 
for a non-vanishing X[^-^ coupling. The region below the dashed curve in (a) is excluded by the L3 
experiment; this limit does not depend on the value of the coupling. Three isolines for the mass of 
the ii are shown as dotted curves. 



couplings were derived. The limits were set within the unconstrained MSSM, the constrained 
MSSM and the minimal Supergravity model. The model dependence of the results was studied 
in detail by performing a scan of the MSSM parameters and was found to be small. 

In the large part of the MSSM parameter space covered by the scan, the existence of squarks 
coupling to an e^d pair with masses up to 260 GeV is excluded at 95% confidence level for a 
strength of the Yukawa coupling equal to a^m- For a 100 times smaller coupling strength squark 
masses below 182 GeV are ruled out. This improves the indirect bounds set by low-energy 
experiments and, in SUSY models where the sfermion and the gaugino sectors are not related, 
extends beyond the reach of other collider experiments. In models where the sfermion masses 
depend on the parameters which determine the supersymmetric gauge sector, the limits extend 
beyond the constraints obtained at the Tevatron collider for intermediate values of tan P and for 
a Yukawa coupling of electromagnetic strength, and are comparable with LEP bounds in part 
of the parameter space. 
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